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SUMMARY 

This work presents experimental data describing the forcing functions and 
structural responses characterizing gas turbine rotor/case system vibration 
due to rotating stall in an axial-flow compressor. Two data sets with funda 
mentally different response characteristics are presented; one is supersynch 
ronous Ind the other subsynchronous. Conventional beam-element rotor dynamics 
analysis is shown to be severely limited in its ability to predict these re- 
sponses. A new analytical approach, which significantly increases structural 
response predictive capability for these phenomena, is briefly discussed. 

INTRODUCTION AND BACKGROUND 

The emphasis on aerothermodynamic performance in the high-performance gas 
turbine engine of today and tomorrow has focused attention on increased per 
formances from its components, one of which is the axial-flow . 

Compressor design trends, aimed toward high efficiencies and/ or high aerody 
namic loadings, have led to machines that are very sensitive to 
clearances between blades/vanes and their respective case/ rotor structures. 

The controlling criterion in the selection of these clearances is often the 
necessity to accommodate rotor and case excursions associated with surge/stall 
instabilities. Closely related to clearance selection is the rotor dynamic 
design of a given machine. The mechanical designer typically has only empiri- 
cal guidelines based on severely limited data with which to design to accommo 
date surge/stall deflections and to simultaneously provide minimum running 
clearances. This paucity of knowledge of structural dynamic response duri^ 
surge/stall instabilities renders the optimization of tip clearances, rotor/ 
case flexibility, and bearing placement a difficult and abstract task. 

Detroit Diesel Allison (DDA) is currently working under Air Force con- 
tract to experimentally document compressor structural response to the rota 

ting stall phenomenon and to improve state-of-the-art ana H^firioncies 

for predicting such response. These efforts address some of the deficiencies 
in kLwledge noted above. This paper presents an abbreviation ^ 

results of this contractual effort with emphasis on the experimental data. 

Rotating stall can be viewed as one or more regions or "cells" of low mo- 
mentum low-pressure fluid, each covering some sector of the flow annulus o 
the compressor. Although rotating stall typically occurs near 
surge line, some machines exhibit it away from surge, even near 
linL The number of cells can vary from one to eight or more. Multiple cells 
to be u„i£o«ly .paced circumfereetUlly. In the folly 
ting stall considered herein, the cells traverse the axial length of the com 

♦This work is associated with U.S. Air Force Contract No. F33615-79-C-2089. 
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pressor, so^^can be regarded as "tubes" of low-pressure fluid that are axially 

£o“« fhese stall calls result in net 

compressor rotor and the case, assuming the cell 
Is the flow passage span. These forces are radial and, if ramp angles are 

tiorof’fh*" components. The stall pattern rotates in'^th^d irec- 

Ind compressor rotor but at partial rotor speed, generally between 25% 

stall 'T a%ingle-c^ll ^"^ti^ 

two-cell stall at ^ subsynchronous excitation mechanism. A 

s^an at Ini A appears as a synchronous force, while a two-cell 

stall at 60% speed appears supersynchronous. 

atPai" ®crge phenomenon, rotating stall can produce quasi- 

steady, periodic forcing functions which can cause steady-state, periodic 
rotor/case response manifested as structural dynamic instabilities. Such re- 
sponses can, based on experience, lead to high amplitude vibration and/or 
radial excursions beyond the design clearances, causing rub and potential 
structural failures. This paper presents a classic set of experLntal data 

bratio^"%i nature of rotating stall-induced rotor/case system vi- 

bration. The associated forcing functions are also presented, and pLrier 
analysis is utilized to gain insight into the forcing function/ response rela- 
tionships. The new analytical approach developed at DDA to calculLe stall- 
induced system response is briefly outlined. ^ 

EXPERIMENTAL DATA AND DISCUSSION 

presented herein were generated at DDA on a large, multistage, 
compressor rig, which is fairly typical of the compressors current- 

ed irFieure ir turbines. This rig is sketch- 

mpnt of fr .Though not shown in this sketch, the bearing/rotor arrange- 
e rig IS similar to that of the parent gas generator, of which this 
i^r^rr component. This rig was instrumented to measure the time-var- 

lant forcing functions and structural responses associated with rotating 

Fillip 2 high-response instrumentation scheme are shown in 

igure 2, which depicts the compressor case unwrapped into a plane. These 
sensors are all mounted on the compressor case in the stationary frame. The 
^namic pressure sensors are high-response transducers flush-mounted in the 
outer gas-path wall. They measure the fluctuating component of endwall pres- 
sure, which IS assumed to characterize the stall cell at both its inner Ld 
outer span. Vibration pickups measure case velocities, while rotor-to-case 
cated^""^ motion is sensed with "whip" pickups, which are proximity probes lo- 

ia?iv HA- Rotor motion can be established by vector- 

laUy adding case motion and rotor-to-case relative motion after transforming 
these measured parameters to compatible units. ^ 

Running of this compressor rig produced several rotating stall and/or 
surge sequences, which were recorded on magnetic tape via the described in- 
strumentation. The phenomenological sequence from which the data presented 

quantified consists of a high-speed surge followed by two-zone and 
then single-zone rotating stall patterns. Figure 3 is a time record of the 
rotor speed excursion showing the onset of surge and the subsequent initiation 
varianf* Patterns. As seen in this figure, the entire sequence of time- 

variant phenomena lasted approximately 5 seconds, during which time the rotor 

fTT rt T" influence o£®the r^«"d 

fwrtW I f I presented in the subsequent discussion were defined over^he 
two time intervals shown in Figure 3. 
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Following surge, s two-rone rotating stall pattern was spontaneously es- 
tablished. traveling (one-per-circuraference spatial propagation rate! a 
rotirTeed Time Laces of the dynamic signals over the defined ana ysis 
interval are presented in Figures 4, 5. and 6. A common timing signal on 
these plots is used to time-relate events. The responses shown in Figures 5 
and 6 show high-amplitude compressor case structural vibration 
quency of the individual stall cell passage from Figure 4, i.e. at two times 
Ihe tL-zone rotating stall rotational frequency. The whip signals ^ 

6 also indicate rotor-to-case relative motion at the same frequency. h 
responses and the associated forcing functions thus appear to be supersynch 

ronous. 

The two-zone rotating stall pattern evolved into a si^le-zone 
(Figure 3) with a propagation rate of 60% rotor speed. Figures 7, 8 and 9 
present the traces of dynamic pressure, case vibration, and whip for this 
stall occurrence. Again, high amplitude structural response is observed. 

Both thrcompressor Ls. respou.. and the whip aignals visually tndtcate major 
coLponLta at the subaynchronoua paaai.« frequency “'.‘''e Ftgure 7 orcr^ 
funLions. Clearly evident in some of the response signals are additional 

higher frequency components* 

Insight into the nature and relationships of the forci^ ^T^^ThrFoi.r- 
structural responses characterizing these stall phenomena is affor ® ^ 

iL analysis. The two-zone stall data of Figures 4 through 6 were digitized 
over the common interval shown on Figure 4, and the “"^le-zone data o g 
ures 7 through 9 were likewise digitized over the period shown on Figure 7. 
Digitized data were harmonically analyzed to produce components for each sig 
nal according to the cosine Fourier form. 


n=0 


Um cos 


( m ft t - md+ a^j) 


where : 

Uju = amplitude of mth harmonic (harmonic coefficient) 
m = harmonic index 
ft = fundamental frequency 
d = circumferential location of sensor 
= phase of mth harmonic 
t = time 


This procedure removes phase differences arising from the different circumfer 
ential locations of the various sensors. Note that in Figures ^ 
representing fully developed stall patterns, all these dynamic signals, co 
siLred to Lpresent transient phenomena in the sense of an operating steady- 
state turboraachine, are relatively periodic functions. 

The Fourier content of the pressure traces is presented in Tables I and II 
for the two-zone and one-zone phenomena respectively. 

circumference spatial cycles were digitized (Figures 4 and 7), the first five 

spatial harmonics of the waveforms correspond '^®”P°’^Lgurrr 

ics 8, 16, 24, 32, and 40, as listed. A comparison of Table 1 

for two-zone stall confirms the predominance of second harmonic con en 
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circumference) and indi- 

lall\t 11 T harmonic. Comparison of Table II with Figure 7 for one- 

zone stall shows a predominant first harmonic content with significant higher 
harmonic content, notably second harmonic* 

Table III presents the Fourier content of the measured responses to the 

J are seen to be axially toward the rear of the compressor, whereas the 
highest amplitude case responses from Table III are seen to be in the compres- 
sor inlet area, suggesting a proximity to a resonance in the inlet. Compari- 
on of these two tables does clearly show, however, direct structural response 
to rotating stall. The presence and predominance of second and higher haLon- 
ic case response show the importance of localized distortion and ••Sging"™" 

ftrsrh^llL r behavior, correspfnd?^ to 

uto^to response components, is a relatively minor con^^ib- 

scrip^Lns Figure 10 further clarifies these de- 

zone^stall^'^ Kramn ? of the measured responses to single- 
Wi?h ! amplitude predominance of first harmonic is apparent, al^g 

with a strong second harmonic contribution, especially in the i^et. This * 

r.r »hiP aignal. The reUt'.; pUjiier 

contents these response signals are comparable to those of the forcing 

functions (Table II), again suggesting the approach of comparing forces and 

ITITIV '’‘T";' PrLary^eaL “epoLe 

be beam like, unlike that observed in the two-zone phenomenon. 

As shown in Figure 2. axially-redundant signals were obtained for several 
p rameters. Examination of the Fourier content of these signals in Tables III 

sor in botrfLsrand^^’'^"''H^i^^ asymmetric responses throughout this compres- 
caL signals ?L - ®-f ° components for both case and rotor-to- 

M,- * he Significance of these observations is that the responses of 

this axisymmetric turbomachine to these stall-induced forcing functions are 
indiV axisymmetric. The asymmetric case response in the^irst harmonic 

indicates asymmetric mounting of the case to ground for the beam-like motion 
Asymmetric case response in the second and higher harmonics is due to asy^- 

metry in the case itself, caused by holes, split lines, bleed manifo^r Ld 
other deviations from pure axisymmetry. niroios, and 

of *^bese experimental data is that stall-induced response 

of the compressor case, an important contributor to system dynamic response 
bv ^°bor-to-case clearance response, cannot always be adequately described 
by first harmonic, beam-like motion. Furthermore, significan? c ircumf^eet i^t 
asymmetry is present in the different modes of case response. Obvious from 
these observations is the fact that the beam-element analytical aepr^Ich wTde- 
If th^ analyses is generally inadequate for modeling either 

two-zLf°8tall’'ir8een described herein. Since the supersyUfhronous 

o zone stall is seen to have only minor components of forcing functions and 

respon... that ca„ be described by the first spatial hanaoniirthrbeam-e la- 
ment approach is inherently grossly inadeqnate. Although forcing functions 
an reaponsea of the subsynchronous single-zone stall have predominant first 

hara»nic content the higher harmonic components are not negligible bmta"e 
Ignored in beam-element analysis. Si-^Kioie duc are 
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NEW ANALYTICAL APPROACH 


DDA has developed a capability for predicting rotor/case 
spouse to dynamic forcing functions. This analysis has ^een developed specif 
icallv for the rotating stall phenomenon in compressors. The analysis in 

iZ abllUy to pradict localised motion, of the at.tio atructura (t.a., 
:i::J“,.or :La) /including ■•agging" modaa. Tha approach 

beam-element rotor analysis along with an axisymmetric finite element case 
analysis. Egging and localized motions are calculated for the case alone 
^^r^sLed to he decoupled from the rotor. Beam-like motion is calculated 
for the coupled rotor/case system using finite-element case modes modal^ 
coupled to the conventional beam-element rotor model. 

axisymmetry are recognized: (1) in the beam-element calculation, through 

’ lowance for asymmetric support modeling, and (2) in the finite-e emen 
.rir.t.;"hrough availability of a diacoutinuoua .lament for modaUng atruta. 
vanes holes in the case, bleed manifolds, and other variances from pure axi 
sy^e^rj! This techniqu; has substantially addressed the shortc^i^s of con- 
vLtional beam-element rotor dynamics analysis in the mode li^ of s^all in 

duLd response and has produced reasonable analytical simulations of both the 

sub synchronous and supersynchronous preceding experimen 
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TABLE I. - FOURIER DECOMPOSITON OF DYNAMIC PRESSURES FROM TWO-ZONE STALL 



Temporal harmonic * 8 


Amplitude 

Phase 

0.214 3 
“103.1 

0.0346 

-75.9 

0.0613 0.0670 0.0562 
97.0 -144.7 135.5 

0.0475 

-149.5 

0.1012 

-174.5 

0.0924 

60.7 

0.2463 

77.7 




Temporal harmonic - 16 





Amplitude 

Phase 

0.1050 

-1.5 

0.1597 

-122.6 

0.1961 0.9671 1.784 

78.1 -160.2 -176.8 

2.711 

-175.6 

2.622 

159.5 

3.878 

136.4 

2.340 

106.1 




Temporal harmonic ■ 24 





Amplitude 

Phase 

0.0498 

47.3 

0.0334 

18.4 

0.0260 0.0513 0.0620 
115.3 -74.0 140.6 

0.0561 

153.8 

0.0688 

166.4 

0.0974 

88.4 

0.1809 

61.6 




Temporal harmonic - 32 





Amplitude 

Phase 

0.1022 

-22.9 

0.1339 

-24.8 

0.102 3 0.1286 0. 3121 
-28.0 -164.8 -171.1 

0.5181 

-154.2 

0.4338 

141.5 

0.8993 

102.6 

0.1992 

117.0 




Temporal harmonic “ 40 





Amplitude 

Phase 

0.0380 

0.0 

0.0026 

-180.0 

0.0001 0,0311 0.0495 

-180.0 0.0 -180,0 

0.0168 

0.0 

0,0526 

0.0 

0.0914 

0.0 

0.0211 

0.0 


Amplitude in pal. 
Phase In degrees . 


TABLE II. - FOURIER DECOMPOSITION OF DYNAMIC PRESSURES FROM SINGLE-ZONE STALL 



Amplitude 

Phase 

0.2323 

17.3 

0.2519 

149.3 

0.3808 

176.0 

1.041 1.576 

-157.3 -172.7 

2.222 
-17 3.5 

2.122 

177.1 

2.139 

155.6 

0.7372 

128.1 




Temporal harmonic ■ 16 





Amplitude 

Phase 

0.1429 

-7.6 

0.0926 

38.3 

0.1158 

-108.3 

0.4098 0.7246 
-134.4 -147 .7 

1.111 

-155.3 

0.9731 

179.7 

0.9519 

166.4 

0.4192 

108.8 




Temporal harmonic ■ 24 





Amplitude 

Phase 

0.1128 

-26.2 

0.0888 

7.0 

0.1165 

^7.2 

0.1704 0.3000 

-146.3 -145.2 

0.4062 

-150.8 

0.3069 

-164.6 

0.2796 

-127.8 

0.004 3 
169.1 




Temporal harmonic “ 32 





Amplitude 

Phase 

0.0399 

-32.8 

0.0845 

32.1 

0.0883 

-32.6 

0.0797 0.1446 
-124.5 -143.3 

0.1212 

176.6 

0.0718 

-40.8 

0.1437 

-30.1 

0.1235 

127.2 




Temporal harmonic • 40 





Amplitude 

Phase 

0.0085 

0.0 

0,0129 

0.0 

0.0245 

0,0 

0.0408 0.0491 
-180.0 -180.0 

0.0655 

-180.0 

0.0659 

0.0 

0.1723 

0.0 

0.0220 

0.0 


Amplitude in psi. 
Phase in degrees. 
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FOURIER COMPONENTS OF EXPERIMENTAL RESPONSE FOR TWO-ZONE ROTATING STALL 
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TE-6657 

Figure 1. - Compressor rig on test stand. 
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Figure 2. - H i gh- response instrumentation scheme. 
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Rotor speed— rpm 



Figure 3. - Sequence of phenomena following surge. 
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Fifure 4. - Pressure 4«ta d^rlnf tw«-zM« stall 



Figure 5. - Vibration induced by two-zone stall 
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Figure 7. - Pressure data during single-zone stall. 





Figure 8. - Vibration induced by single-zone stall 
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Figure 9. - Vibration induced by single-zone stall. 
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Figure 10. 


Case motion 


descriptions. 
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